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Abstract 

We investigate the Casimir force acting between real metals at nonzero tempera- 
ture. It is shown that the zero-frequency term of Lifshitz formula has interpretation 
problem in the case of real metal described by Drude model. It happens because the 
scattering theory underlying Lifshitz formula is not well formulated when the dielectric 
permittivity takes account of dissipation. To give the zeroth term of Lifshitz formula 
the definite meaning different prescriptions were used recently by different authors with 
diversed results. These results are shown to be improper and in disagreement with ex- 
periment and the general physical requirements. We propose the new prescription 
which is a generalization of Schwinger, DeRaad and Milton recipe formulated earlier 
for ideal metals. On this base the detailed numerical and analytical computations of 
the temperature Casimir force are performed in configuration of two plane plates and a 
spherical lens (sphere) above a plate. The corrections due to nonzero temperature and 
finite conductivity found in the paper are in agreement with the limiting case of perfect 
metal and fit all experimental and theoretical requirements. Among other facts, the 
previous results obtained in frames of plasma model are confirmed. It appears they are 
the limiting case of Drude model computations when the relaxation parameter goes to 
zero. The comparison with the Casimir force acting between dielectric test bodies is 
made. 
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I INTRODUCTION 



The H.B.G. Casimir effect [y predicted more tfian fifty years ago is one 
of the most interesting manifestations of zero-point vacuum oscillations of 
quantized fields. The Casimir effect implies that there is some force acting 
between two uncharged bodies closely spaced in the vacuum. This effect 
is purely of quantum origins. There is no such force in classical physics. 
Unique to the Casimir force is its strong dependence on shape, switching 
from attractive to repulsive as a function of the size, geometry and topology 
of the boundary. The force results from the alteration by the boundaries of 
the zero-point electromagnetic energy that pervades all of space as predicted 
by quantum field theory. Alternatively, the Casimir force can be described as 
the retarded electromagnetic interaction of atomic and molecular dipoles and 
was extended to forces between macroscopic dielectric bodies characterized 
by some dielectric constant [)2|. 

In recent years the Casimir effect has attracted much attention because 
of numerous applications in quantum field theory, atomic physics, condensed 
matter physics, gravitation and cosmology, and mathematical physics (see 
monographs [0, 0, 0] and references therein). New precision experiments 
have been performed on measuring the Casimir force between metallic sur- 
faces I?], g, 0, [TO], |ri], [T2j]. In Refs. |T3|, [T^, [T5j] some promising apphcations 
of the Casimir effect were proposed for diagnostic in thin films and in nano 
electro-mechanical systems. Given the above reasons it is very important to 
understand the Casimir force between real materials including the effect of 
such influential factors as surface roughness, finite conductivity, and nonzero 
temperature. 

Finite conductivity corrections to the Casimir force have long been in- 
vestigated. They were calculated using plasma model |TB|, [T7], jlBj up to the 
first order in the relative penetration depth of the zero-point electromagnetic 
oscillations into the metal. In [|TD| more exact results up to the second order 
were obtained, and in — up to the fourth order. In Refs. [[2T| , |22|] the finite 



conductivity corrections to the Casimir force were computed using tabulated 
optical data for the frequency-dependent complex refractive index. For all 
separations between the test bodies larger than the effective plasma wave- 
length of the metal under study the results of and [jTT], are shown to 
be in good agreement. The effects of surface roughness in combination with 
finite conductivity were investigated in detail in for the configuration of 
a sphere above a plate used in experiments 0, [TD], [n|. Roughness con- 
tributions to the Casimir effect (including origination of a lateral force) has 
been treated recently in [2^, |27| . 

The action of nonzero temperature on the Casimir force between dielectric 
semispaces was taken into account in the Lifshitz theory []2|, |rB|. For perfect 
metals at nonzero temperature the Casimir force was calculated in []2^, 
within the limits of quantum field theory in terms of the free energy density 
of vacuum. There were apparent differences between the results of [@, |T^ and 
which were resolved in [|T8j. As shown by Schwinger, DeRaad and 



Milton [|TBl, to obtain the case of perfect conductor from the Lifshitz the- 
ory one must take the limit of infinite dielectric permittivity before putting 
the frequency equal to zero in the temperature sum. The results of []2|, [TBj 
adapted for the case of perfect metal and of ^ are then in agreement. 
The temperature corrections to the Casimir force turned out to be negligible 
in the experiments 0, [T^, [TT| where the measurements were performed in 
the separation range a < 1 fim. However, at a > 1 fim, as in [|7|, the tem- 
perature corrections make large contributions to the zero-temperature force 
between perfect conductors (e.g., for a = 5 fim the temperature correction in 
configuration of a spherical lens above a plate exceeds the zero-temperature 
force [pOf ). 



The increased accuracy of the Casimir force measurements invites further 
investigation of the temperature corrections in case of real metals. Although 
from a conceptual point of view the Lifshitz theory provides a way of obtain- 
ing all the required results, the problem here is worse than it was with the 



case of perfect metals. In it was suggested to use the plasma model in or- 
der to describe the dielectric permittivity along the imaginary frequency axis 
in the Lifshitz formula for the Casimir force at nonzero temperature (note 
that in [jTB], [T7|, |T^, [T^, ^ the plasma model was apphed for this purpose at 
zero temperature only). In Refs. ^ the detailed calculations of the tem- 
perature Casimir force were performed in framework of the Lifshitz theory 
and the plasma model. It was shown that the temperature corrections are 
negligible at small separations where the finite conductivity corrections are 
very important. By contrast, at large separations finite conductivity correc- 
tions can be ignored, whereas the temperature corrections play an important 
role. In [|32| , p3| the transition region between these two asymptotic regimes 
was also investigated where the combined effect of nonzero temperature and 
finite conductivity is important and should be taken into account. 

It is a common knowledge that at small frequencies the dielectric per- 
mittivity £ is proportional to cj^^ as is given by the Drude model. Because 
of this, the Drude model is favoured over the plasma model (which implies 
£ ~ cj^^) when calculating the Casimir force at nonzero temperature. The 
first attempts to calculate the nonzero temperature Casimir force between 
real metals based on the Lifshitz formula and Drude model were undertaken 
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p3| , p^ , |35| , |36| , pTf . They have led to distinct and unexpected results. 



It was found that the value of the transverse reflection coefficient r2 of the 
electromagnetic oscillations at zero frequency becomes indefinite when one 
describes the boundary made of real metal using the Drude model. In |35j 
the value r2 = has been adopted. This assumption leads to a nonphysical 
conclusion that the asymptotic Casimir force at high temperature in the case 
of real metals is two times smaller than for the case of ideal metal (without 
regard to the particular value of the conductivity of real metal). Also, as a 
result of assumption, made in there arise large negative temperature 

corrections at small separations which are linear in temperature. These cor- 
rections are not only unacceptable from the theoretical point of view but are 



in a conflict with experimental data 



The other authors assumed the value \r2\ = 1 at zero frequency 

using the prescription of [|TB[ formulated for perfect metal and well known 
relation by Hagen and Rubens which is valid for real photons only 
This assumption also leads to nonphysical conclusions, i.e. to linear (although 
positive) temperature corrections at small separations and to the absence of 
any finite conductivity corrections to the Casimir force for real metals starting 
from the moderate separations of several micrometers regardless of metal 
quality (note that the same assumption was accepted in the latest version of 



The situation was clarified in where the discontinuity of the trans- 
verse reflection coefficient as a function of frequency and photon momentum 
was demonstrated in the case of real metal described by the Drude model. 
According to to clear away the ambiguity in the zero frequency term 
of Lifshitz formula arising due to this discontinuity it is necessary to use an 
alternative representation for it [jTB], ^ . This representation gives the possi- 
bility to redefine the zero-frequency term of Lifshitz formula in order to assign 
it the definite meaning for real metals in accordance with the usual physi- 
cal requirements. In however, the values of the Casimir force including 
nonzero temperature and finite conductivity in the framework of the Drude 
model were not computed. Thus, up to now, there is no plausible qualitative 
information on the Casimir effect between real metals including the tem- 
perature corrections. The need for such investigation is apparent when the 
experimental and technological applications of the Casimir force mentioned 
above are considered. It is also important that recently, the Casimir effect 
has been used to obtain stronger constraints on the constants of long-range 
interactions (including corrections to the Newtonian gravitational law) pre- 
dicted by the unified field theories, supersymmetry, supergravity and string 
theory [j3D], [53], The reliable theoretical values of the Casimir 



force at nonzero temperature between real metals and the extent of their 



agreement with experiment are of particular interest in order to obtain the 
strongest constraints. 

In the present paper we propose a method which allows one to attribute a 
definite value to the term of Lifshitz formula at zero frequency for real metals 
as described by the Drude model. This method avoids the above mentioned 
contradictions and solves the problem in a physically consistent way. The 
detailed computations of the Casimir force between real metals at nonzero 
temperature are performed for the configuration of two plane parallel plates 
and a sphere (lens) above a plate in a wide separation range from 0.1 fim to 
10 fim. It is shown that at small separations (low temperatures) the tem- 
perature corrections are small irrespective of whether the Drude or plasma 
model is used. In particular, no corrections to the force arise which would 
be linear in temperature, such as in At large separations 

(high temperatures) there is some difference between the finite conductivity 
corrections to the temperature Casimir force obtained with the plasma and 
Drude models, although the corrections themselves are rather small. With 
decreasing relaxation parameter the results from both models coincide. Our 
results at large separations join smoothly with the increase of conductivity 
with the asymptotic values obtained for the perfect metal (this is not the case 
|35|j where the asymptotic force for a real metal is two times smaller 



m 



than for the perfect one). Also the non-physical results of PT| , according 
to which at large separations the finite conductivity corrections are absent in 
the case of real metals, are shown to be in error. Below the nonzero temper- 
ature Casimir force between the dielectric bodies is also computed and the 
distinctions between the cases of metallic and dielectric bodies are discussed. 
The obtained results are in agreement with the performed experiments. They 
establish the theoretical basis for new precise experiments on measuring the 
Casimir force. 

The paper is organized as follows. In Sec. II the general formalism for 
the Casimir force between real metals at nonzero temperature is presented 



for the configuration of two plane parallel plates. Here special attention is 
paid to the indefinite character and discontinuity of the zeroth term of the 
Lifshitz formula in case of the Drude model. The representation of this term 
is given which is in accordance with the general physical requirements. In 
Sec. Ill the same is done for the configuration of a sphere (lens) above a 
plate. Sec. IV is devoted to the numerical and analytical computations of the 
temperature Casimir force for the configuration of two plane parallel plates. 
Both the asymptotics of low and high temperatures are considered and also 
a transition region between them. Sec. V contains the results of analogous 
computations for the configuration of a sphere (lens) above a plate. In Sec. VI 
the case of dielectrical test bodies is considered and the Casimir force at 
nonzero temperature is found. In Sec. VII reader will find conclusions and 
discussion. Appendices I and II contain some details of the mathematical 
calculations. 

II TEMPERATURE CASIMIR FORCE BETWEEN 
TWO PLANE PARALLEL PLATES 



The original Lifshitz derivation [)2|, |T^ of the temperature Casimir force be- 



tween two semispaces was based on the assumption that the dielectric ma- 
terials can be considered as continuous media characterized by randomly 
fiuctuating sources. The correlation function of these sources, situated at 
different points, is proportional to the ^-function of the radius-vector joining 
these points. The force per unit area acting upon one of the semispaces was 
calculated as the fiux of incoming momentum into the semispace through 
the boundary plane. This fiux is given by the appropriate component of the 
stress tensor (zz-component if xy is the boundary plane). Usual boundary 
conditions on the boundary surfaces between different media were imposed 
on the temperature Green's functions. To exclude the divergences, the values 
of all the Green's functions in vacuum were subtracted from their values in 



dielectric media. The other derivation of Lifshitz formula is based on the 
solution of Maxwell's equations with appropriate boundary conditions on the 
surfaces separating different media In this manner the allowed 

surface modes and a harmonic oscillator free energy of each mode can be de- 
termined. The total renormalized free energy is obtained by the summation 
over all modes and subtraction of the free energy of vacuum. The nonzero 
temperature Casimir force is finally given by the negative derivative of the 
renormalized free energy with respect to the distance between semispaces (at 
zero temperature this procedure is presented in detail in |17]]). 

The modern derivation of Lifshitz formula [j55| is based on the tempera- 



ture field theory in Matsubara formulation. In case of a static Casimir effect 
the system is in thermal equilibrium. To describe it at nonzero temperature 
one should take the Euclidean version of the field theory with a field periodic 
in the time variable within the interval (5 = h/{kBT), where T is tempera- 
ture, kB is the Boltzmann constant. To find the free energy in case of two 
semispaces the scattering problem on the z-axis perpendicular to the bound- 
ary planes is considered. An electromagnetic wave coming from the left in 
the dielectric is scattered on the air gap and there exist a transmitted and 
a reflected wave. Finally the free energy and the Casimir force is expressed 
in terms of the scattering coefficient on the imaginary axis. Calculating this 
coefficient for the problem of two semispaces (notated by an index ss) with a 
frequency dependent dielectric permittivity £{io) separated by a gap of width 
a one finally obtains the Casimir force in the form 



Fss{a) = — ^ J k±dk±qi\ ri^{^i,k±)e 
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where ri 2 are the reflection coefficients with parallel (perpendicular) polar- 
ization respectively given by 



g(^6)g/ + ki 



qi + h 
qi - h 



(2) 



Here is the momentum component lying in the boundary planes, k±_ = 
io = lE,., and the following notations are used 



(li = 



+ kl, ki = 



27r/ 



Introducing a new variable p according to 
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(3) 



(4) 



we rewrite Eq. (|I]) in the original Lifshitz form [j^, |rB|, 



FJa) = -'^t'(:}lp'dv 
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where 



Km-p^ 



Km 



1 



1/2 



(5) 
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and the prime near the summation sign means that the zeroth term is taken 
with the coefficient 1/2. 

Note that the representation of Eq. (|^) for the nonzero temperature Casimir 
force has a disadvantage as the / = term in it is the product of zero by a 
divergent integral. This is usually eliminated by introducing the variable 
y instead of p 



y = 



6 



P. 



6 = 2a^, 
c 



where is the dimensionless frequency. In terms of new variables Eq. 



(7) 

is 
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r2(6,2/) = 
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\e - i)g + y ' 
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(9) 



y + ^{e-l)if + y^ 



£ = em = ^ I 



Both changes of variables and (|7|) are, however, singular at / = 0. In fact 
(@) can be obtained from (|I]) by the regular change of variable 

4a'kl = y'- if. (10) 

Because of this, the equivalent representations of Lifshitz formula ([!]) and (g) 
are preferred as compared to (§). 

To calculate the Casimir force at nonzero temperature between real metals 
one should substitute the appropriate values e{i^i) into Eq. ([TD or (g). In some 
frequency range £{i^i) can be found by the use of optical tabulated data for 
the complex refractive index for the metal under consideration tabulated in, 
e.g., [j5r| (this was done in []2T|, ^). But in any case the optical data should 
be extrapolated outside the region where they are available in the tables to 
smaller and larger frequencies. This can be performed by the use of plasma 
model function 

= 1 + p (11) 

or more exact Drude model one 

where tUp is the plasma frequency and 7 is the relaxation frequency. For some 
metals, the model dielectric functions with appropriate values of cOp and 7 can 
be reliably used throughout the whole spectrum. Calculations of this type 
were performed in [j^, [33], 0, [36], §7^ with different results. As shown 
below, the reason for this is the discontinuity of the transverse reflection 
coefficient r2 at zero frequency with respect to relaxation parameter which 
must be made continuous in a physically reasonable way in order to describe 
real metals by the Lifshitz formula. 

If we put = in Eqs. (g), (|^) and use the plasma model dielectric 
function (p) the result is ^ 

2 



r?(0,A:^) = l, rl{0,k^] 



^k^-Jkl + p 



(13) 



In the case of the Drude model dielectric function ([T2]) it holds 



= 0, (14) 



which leads to 135 



rf(0,A:^) = l, ri{0,k^) = (15) 

for any k± ^ 0. Note that (|r5j) is valid for arbitrary small 7 and arbitrary 
large cjp. Thus the second equation from ( [151 ) is in contradiction with the 
limiting case of ideal metal r^O, ky) = 1 which follows from ( [C3|) in the limit 
^ 00. The results obtained using (fr51) are in conflict with the known 
results for ideal metal [TBI, What is more important, the second equation 
of ([15|) does not approach the second equation of (|I3j) when the relaxation 
frequency 7 goes to zero. In the limit 7^0 one still has r|(0, /c^) = 0, not 
(P^), although XmiEB = in accordance with Eqs. ([TT]), ([T^). In Sees. IV, V 
the discontinuity of the zero frequency term of the Lifshitz formula with re- 
spect to the relaxation parameter will be discussed in detail when comparing 
the continuous modification of this term as suggested below. Exactly the 
same results, as in ([13|), ([T^), are obtained if one uses the variables (^, y) 



instead of (^, k\}) and the representation (^ of the Lifshitz formula. 

If the original representation (|^) of the Lifshitz formula is exploited and 
written in terms of the variables (^, 'p) the situation changes drastically. Here 
one immediately arrives at 

T\{^,V) = rl{^.V) = ^ (16) 

for all j3 7^ 00 irrespective of whether the plasma or Drude model is used for 
the dielectric permittivity on the imaginary axis. The reason for the distinct 
value of r2(0,j9) is the singular character of the change of variables (|5) when 
1^0 = 0. This change relates the single value of /c^ = to all the finite values 
of ]7, and all the values of /c^ 7^ to a single point p = 00. The values ([TC| ) for 
both reflection coefficients at zero frequency were postulated at [BB, p7j and 



used in all the numerical computations even after one more singular change 



of variables (|7|) was performed in order to calculate the divergent integral in 
Eq. (|5]). We show below, that this postulate, though it works good for the 
perfect metal, is not justified for real metals of finite conductivity. 

In actual truth, the transverse reflection coefficient r2 from Eq. (0) or 
(0) has a discontinuity as a function of two continuous variables and y 
for instance) in the case of the Drude dielectric function [j33j. If one puts 
^ = in (1^) from the very beginning the result is r2(0,?/) = in agreement 
with (p^). If, however, one approaches the point = 0, = 0) along the 
direction ^ = ky in the (^,?/)-plane then r^H^y) 1- The distinguishing 
feature of the Drude model is that for an arbitrarily small y there exists k 
such that r2{ky^y) takes any value in between zero and unity. By contrast, 
in the case of plasma model both reflection coefficients are continuous, and 
r2{ky,y) does not depend on k but is determined only by the value of y. 
Of even greater concern is that the transverse reflection coefficient at zero 
frequency is discontinuous with respect to the relaxation parameter 7. 

As is seen from the above discussion, there is a serious unresolved issue 
concerning the value of the zero frequency term of the Lifshitz formula for 
real metals. In fact the scattering problem, which forms the basis of the 
Lifshitz formula, is meaningful only for nondissipative media (this is the case 
for plasma model or for dielectric materials). The case £ = 00, as it is for 
metals at zero frequency, is especially complicated when the Drude model is 
used taking account of dissipation. The latter leads to the violation of the 
unitarity condition and thereby the zeroth term of Lifshitz formula becomes 
indeflnite and must be redeflned. For ideal metals of inflnite conductivity 
this issue was resolved by the Schwinger, DeRaad, and Milton prescription 
[ [TBI demanding that the limit e ^ 00 should be taken before setting ^ = 
(which is equivalent to the use of Eq. (|TBO). Let us flnd out, in what way the 
problem of zero-frequency term of the Lifshitz formula can be solved for real 
Drude metals in a physically satisfactory way. 

For this purpose let us use the representation of the Lifshitz formula in 



terms of continuous frequency variable. Such a representation was suggested 
2B|] and used in [jT^, ^ for the other purposes. According to the Poisson 
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summation formula if c{a) is the Fourier transform of a function b{x) 

-| oo 

c{a) = 1^ J h{x)e-''^'' dx (17) 



27r 

then it follows 



-oo 



00 



E h{l) = 27r E c(27r/). (18) 

l=—oo l=—oo 

We apply this formula to Eq. (g) using the identification 



\1\t 

where = rl and 

fss{lr,y) = fi^{lr,y) + fi^{lT,y), (20) 
/i?(/r,^) = (rrV-l)-\ z = 1, 2 



are the even functions of /. 

Then the quantity Css{a) from Eq. ( \n\ ) is given by 

oo oo 



IGTT^a^ 



J dx cos ax y y'^ dy fss{xT,y). {21] 



Using Eqs. (§), (|TB|), (^ one finally obtains the representation of the 
Lifshitz formula 



Fss{a)= E hs{l) (22) 

l=—oo 

-^^t^dl COS (lilf) J yUyUly), 



where the continuous frequency variable ^ = rx, and kBT^ff = he/ {2a). 
Note that in the representation ([22| ) the / = term gives the force at zero 
temperature. It is useful also to change the order of integration in Eq. ( |22| ) 



Fss{a) = -jTr^.T. j y'dy / d^ cos {l^t) fU^,y), (23) 



oo 
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where t = T^ff/T. 

Let us isolate the zero-frequency term of the usual Lifshitz formula (§) 
with discrete frequencies in representation (|23| ). In this a way it will be 
expressed in terms of the integrals with respect to continuous variables ^, y. 
For this purpose we write out separately the terms of Eq. (^B)) with / = 0, 
integrate all the other terms by parts with respect to ^, change the order of 
summation and integrations, and use the formula 

sm{lzt) 1 
~ 2 



1=1 



I 



TT + 27tA (^] 



\27tJ 



tz 



(24) 



where A[z) is the integer portion of z. The result of these transformations is 



Fss{a) = 



he 



327r2a4 



y^dy / dUss{L y) 
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y^ dyfss{y,y) 



TT + 27rA 
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TT + 27CA I — 
27T 



Here the first term is the zero-temperature force. The second and the third 
contributions in the right-hand side of Eq. ([25| ) can be transformed using the 
definition of the function A{z) and representation of fss in terms of the sum of 
parallel and transverse modes in accordance with ([2D| ). Taking into account 
that f^l^ and its derivative are continuous functions (as distinct from /j^^) 
one obtains finally after the cancellation of zero-temperature contribution 
(see Appendix I for details) 



a] 



E / y dyfssiiuy)- 



(26) 
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In this equation all the terms with / > 1 coincide with those in Eq. (|^) and 
the term with / = is given by 



a] 



167ra^ 
oo y 
y'^dy I di 





y'dy [fil\o.y) + fS\y:y) 
dfSKly) 



(27) 



The obtained representation for the zero-frequency term of the Lifshitz 
formula is well suited for solving the problem formulated above. In terms of 
dimensionless variables y) the plasma and Drude dielectric functions take 
the form 



~ 2 ~ 2 

r iKi + t) 



(28) 



where cD^ = 2aijjp/c^ 7 = 2a7/c. Evidently, no discontinuity problem arises 
in the single integral with respect to y in the right-hand side of Eq. ( |27|) . It 
is well defined for both dielectric functions of Eq. (|28|) and for both polariza- 
tions, and contains the limiting cases of an ideal metal and zero relaxation 
parameter. What this means is that the single integral term in ( P7|) is an- 
alytic with respect to 7, i.e. its values calculated with the Drude model 
approach the values calculated with the plasma model when 7^0 (see also 
the numerical computations of Sec. IV). 

A completely different type of situation occurs in the double integral from 
the right-hand side of Eq. (|27|). In fact, by the use of (|9D, ( PD| ) one obtains 



= 2- 



where 
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(31) 



'{8 - i)e + 2/2 + ^ 

Substitution of e = £p from Eq. (^) into (p9|), (|30|) leads to 

dr2{ly) _ dfS\ly) _^ 

di di 

so that in the case of the plasma model the double integral in ( |27|) vanishes. 
If, however, e = is substituted into (|2^), ( |3I]|) one discovers a discontinuity 
in 7. Actually, as it follows from (|2^), ( PD|) 

lim ^-f^SK^^y) ^ Q /32) 
7-0 



for any ^ > and y 0. On the other hand, 
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.y - y^p + y' 

To conclude, the presence of the double integral in the right-hand side of 
Eq. (^7|) renders it unsuitable for real metals described by the Drude model. 
To remedy the situation let us notice that under the increase of e it holds 

dfSKly) 



1 



^0 with £ oo. 



(34) 



Consequently, the discontinuous term can be deleted by the use of Schwinger, 
DeRaad and Milton-type prescription [jT^ which is that the limit of dielectric 
permittivity to infinity must be performed before setting frequency to zero 
when considering the zero-frequency term of the Lifshitz formula. For ideal 
metals the use of such a prescription is necessary in order to achieve agree- 
ment between the Lifshitz formula for dielectrics and field theoretical results 
obtained for a perfect metal, i.e. for boundaries with the Dirichlet boundary 
condition. As it is advocated here, for real metals the analogous prescription 
must be used in a more restricted way to eliminate the discontinuity in the 
zero-frequency term. 

With the above prescription the Casimir force acting between two plane 
parallel plates at nonzero temperature is given by 



FJa) = 



167ra^ 
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In Sec. IV the detailed calculations of the Casimir force using Eq. (^ are 
performed. The obtained results are compared with those obtained by other 
authors. 



Ill TEMPERATURE CASIMIR FORCE IN CONFIG- 
URATION OF A SPHERICAL LENS (SPHERE) 
ABOVE A PLATE 



In this section the most important results of Sec. II are adapted for a lens 
(sphere) above a semispace (plate) which is the configuration most often used 
in experiments [)7|, 0, [T^? Ol • Although there is no fundamental derivation 
of the Casimir force acting between a lens and a plate the proximity force 
theorem []52| gives the possibility to express it in terms of the Casimir free 
energy density for the configuration of two plane parallel plates. This latter 
is given by [|32|, [33|, p| 



Essia) = 



47r 



E 



} 



(36) 



dk^{\xi [l-r?(^;,/c^)e"2"'^' 

+ ln [l-rl{iuk^)e'^^^[ 

where the reflection coefficients and other notations are introduced in Eqs. (0), 
(0). Note that the Casimir force between plates from Eq. (|T]) can be obtained 
as Fss = —dEss/da. 

According to the proximity force theorem the Casimir force acting 
between a lens (sphere) of radius R and a plate is 

F,i{a) = 2TiRE.,s{a) 
keTR 



E jk^dk^[\n[l-rl{^i,k^)e~'""i^ 

l = — 00 Q 

+ ln [l-r|te,fc±)e-''""]} 



(37) 



Index si here stands for semispace-lens. This formula is valid with rather 
high accuracy of about a/R []5^ which is usually a fraction of a percent. 

Introducing the dimensionless variables and y from Eqs. (|7|), ( pXI| ) we 
rewrite Eq. (|371) in the form 

keTR ^ 



Fsiia) = ^ E / 2/0^2/ {in [l - rl{ii, y)e 



+ ln [l-rl{i,y)e-y]}. 



(38) 



This result is in direct analogy with Eq. (g) for two plates. There is no need 
to use singular changes of variables as discussed in the previous section. 

The term of Eq. (^) with / = suffers exactly the same ambiguity as 
the zeroth term of Eq. (g) if one uses the Drude model ([T^) to describe the 
dependence of a dielectric permittivity on frequency. This ambiguity lies 
in the fact that the quantity r^i^^y) has a discontinuity at a point (0,0) as 
a function of two variables and is discontinuous also at the point of zero 
relaxation parameter (see the detailed discussion in Sec. II). Because of this, 
the zeroth term of Eq. when applied to real metals, must be modified 
in the same way as it was done for configuration of two plane plates. For this 
purpose we apply the Poisson summation formula ([T^) and rewrite Eq. ( |35| ) 
in the form analogical to Eq. 



Fslia) = ^ E7^^ y^^) I ydyfsi{U). (39) 

where 

fsi{ly) = f^^\ly) + f^^\ly). (40) 



Evidently, the zeroth term of Eq. (|39j) gives us the Casimir force at zero 
temperature. The terms with / > 1 represent the temperature correction. 
After changing the order of integration in Eq. ( |39| ) one obtains 

%cB °° , °r \ ~ 
Fslia) = ^J: ydy d^cos (l^t) U^,y). (41) 

Using Eq. ([2^) and performing exactly the same transformations as in 
Sec. II and Appendix I, the Casimir force between a plate and a lens takes 
the form 

'T^ P? oo ^ 

Fslia) = Ft'\a) + ^ E / ydyUii. v)- (42) 
Here all terms with / > 1 are exactly the same as in Eq. (§§), whereas the 



zeroth term is given by 
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ydy [fij\0,y) + fS\y,y) 

dfS\ly) 



(43) 



This representation of the zero-frequency term is in direct analogy with 
Eq. ([27| ) for two plates. Like in (P7|), no discontinuity is contained in the 
single integral with respect to y. As to the double integral in (|?3|), it is dis- 
continuous with respect to the relaxation parameter 7 at a point 7 = 0. In 
fact, it follows from Eq. (gO|) 



(44) 



di eV-rU^y) 0^ 

where the derivative of the transverse reflection coefficient was calculated in 
Eq. (|3I]| ). Once again, in the case of plasma model the derivative ( |44|) , and 



also the double integral in 
from (mi) 



vanish. As to the Drude model, it follows 



7-0 di 
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(45) 



(46) 



which means the discontinuity of the double integral in (|43|) as a function of 
7- 

This discontinuity leads to the nonphysical results (see Sees. IV, V). It can 
be deleted by taking the limit of the dielectric permittivity to infinity prior 
to setting the frequency equal to zero in accordance with prescription used 
in Sec. II. Then due to (^), 

dfS\ly) 



1 



di 

and the discontinuous term in (|43| ) disappears 



with £ ^ 00 



(47) 



Using this prescription, one finally obtains the expression for the Casimir 
force acting between a plate and a spherical lens (sphere) at nonzero temper- 
ature 



keTR 
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ydy f^l\o,y) + fS\y,y) 







+2J: ydyf,,i(li,,y) 



(48) 



The calculations of the Casimir force on the base of Eq. (|48|) are contained in 
Sec. V. The obtained results are compared with those obtained without the 
use of the above prescription [jS^, p5| or using another prescription [pO, B7I . 



IV COMPUTATIONS OF THE CASIMIR FORCE BE- 
TWEEN TWO PLATES MADE OF REAL METAL 



We start the numerical computation of the Casimir force with Eq. (^5|) using 
the Drude and plasma dielectric functions (|r2|). For the sake of definiteness, 
we consider the case of aluminum whose plasma frequency and relaxation 
parameter can be chosen as follows [|5D| 



cjp « 12.5eV« 1.9 X lO^^rad/s, (49) 
7 « 0.063 eV « 9.6 x 10^^ rad/s. 

The most descriptive quantity illustrating the temperature correction to the 
Casimir force is 



Li J 



Flia ,To)-Fl ia,0) 
F/.(a,0) 

Here the upper index / near the quantity Fss from Eq. (|351) signifies one or 
another dielectric function (Drude or plasma for instance), and the second 
argument Tq = 300 K, and T = means the temperature at which the Casimir 
force is computed. The relative temperature correction ( |5I]| ) is a function of 
separation a. It incorporates effectively both the case of low temperatures 



(which occurs when Tq <C T^fj at small a) and high temperatures (as Tq 
Tg// at large a). 

The results of numerical computations by Eqs. (|55|), (^Dj) the case f = D 
(Drude model) are presented in Fig. |T] by the solid curve 1. In the same figure 
the sohd curve 2 illustrates the result obtained for dielectric plates (this case is 
discussed in Sec. VI). The dashed curve represents the result obtained by the 
approach of In accordance with this approach, as was discussed in 

Sec. II, one should substitute f^s\0, y) = instead of f^fiy, y) into Eq. (|35|). 

As is seen from Fig. |I] (curve 1), in the case of real metals the relative 
temperature correction is rather small at small separations and monotoni- 
cally increases with the increase of separation remaining positive as it should 
be from general thermodynamical considerations. In contrast to this, the 
temperature correction, computed on the basis of |35], is negative in a 
wide range of separations and turns into zero at a separation a « 6.3 fim. 
For larger distances it becomes positive. At small distances it is rather large 
by the modulus and increases linearly with distance (temperature). This 
behavior is in radical contradiction with both the case of ideal metal and a 
metal described by the plasma model dielectric function considered in detail 
. At the same time the results obtained here are in good agreement 



m 



with both limiting cases of an ideal metal and metal described by the plasma 
model. 

Note that in the original paper not 6t{F) but the correction factor 
was plotted as a function of separation a, i.e. the ratio between the Casimir 
force for real metals computed at Tq = 300 K and the Casimir force for ideal 
metals at zero temperature. According to the results of [j32j in the case of 



plasma model this correction factor is approximately equal to the product 
of the correction factors through finite conductivity (at zero temperature) 
and through nonzero temperature (for the ideal metal). Thus the resulting 
negative temperature corrections, which are unacceptable from the thermo- 
dynamical point of view, were not clearly evident in . 



In Table 1 the values of the relative temperature corrections are given 
at several distances between the plates in the case of Drude model (present 
paper, column 2), plasma model (present paper, column 3), for the ideal 
metal (column 4), computed on the base of in frames of Drude model 
(column 5), and computed on the base of ^ in frames of Drude and 
plasma models (columns 6 and 7 respectively). We remind that in |37 



all computations were performed by Eq. (|5|) where in the zeroth term the 
reflection coefficients (|T6|) were substituted. As is seen from the comparison 
of columns 2-4, in the framework of the Drude model the relative temperature 
correction at smallest separations is three orders larger than with the plasma 
model and four orders larger than for the ideal metal. Note that for the ideal 
metal at small separations the lowest nonzero temperature correction to the 
zero-temperature result is of order {T/Teff)"^ (see below). With increasing of 
separation the difference between the predictions of Drude and plasma models 
decreases approaching the values obtained for the ideal metal (column 4). At 
small separations the modulus of temperature correction calculated on the 
base of (column 5) is several times larger than calculated by us (column 
2) and the correction itself is negative. The results of the calculations on 
the base of ^ (columns 6, 7) are also significantly different from our 
corresponding results of columns 2, 3. According to [j36|, ^ there exist large 
temperature corrections at small separations which are linear in temperature. 

Now let us discuss the finite conductivity corrections to the Casimir force 
at nonzero temperature. They can be characterized by the quantity 

*"*^»»'- f»{a,T„) • 

where, as above, F^g{a, Tq) is the Casimir force between perfect metalfic plates 
at temperature Tq. The results of numerical computations of (5c(F/g) in the 
case of the Drude model (/ = D) are shown in Fig. 0. The solid curve rep- 
resents the finite conductivity correction obtained by our method (Eq. (|351) ) 
at Tq = 300 K. For comparison the corresponding results at zero tempera- 
ture Tq = are given by the short-dashed curve. The long-dashed curve is 



obtained at Tq = 300 K in the framework of the approach of In 
Fig. the results computed in the case of the plasma model (/ = p) are 
shown with the same notations as in Fig. @. As is seen from Figs, g, 0, in the 
Drude model the finite conductivity correction is larger than in the plasma 
model at all temperatures and at all separations. From Fig. it follows that 
at a = 1 fim the finite conductivity correction computed using the plasma 
model at Tq = 300 K is practically the same as at Tq = 0. 

It is notable that the finite conductivity correction computed on the basis 
of 13^, ^ turns into zero at separations larger than 5 fim in both Drude and 
plasma models (see the long-dashed curves in Figs. and |^). In fact, as 
is seen from the solid and short-dashed curves in Figs. 0, with increasing 
temperature the finite conductivity correction decreases. The extent of this 
decrease depends on the model used {D or p) and on the values of plasma 
and relaxation frequencies. In contrast to this, the finite conductivity correc- 
tion computed on the basis of ^ turns into zero at one and the same 



separation not only in different models but also for one and the same model 
with different values of parameters, i.e. regardless of the quality of metal 
under consideration. It is easily seen, that the separation value, at which the 
finite conductivity corrections of ^ vanish, is determined not by the 
conductivity properties of a metal but by the smallness of terms with / > 1 
in the Lifshitz formula which are of order exp(— 27rr/re//). This means that 
the anzats used in [|35i |37|] to redefine the zeroth term of Lifshitz formula (see 
the above discussion in Sec. II) is unjustified. 

It is of interest to consider the dependence of our results on the relaxation 
parameter 7 of Drude model. In Fig. g the relative temperature correction 
computed by Eqs. (^5]), ( |5D| ) is plotted with the decreasing of 7 at a separation 
a = 2 fim, and at a temperature Tq = 300 K (solid curve). In the same figure 
the dashed curve represents the result obtained using the plasma model. It 
is clearly seen that with the decrease of 7 by 3 orders of magnitude the 
results computed in the Drude model join smoothly with the results of the 



plasma model as it should be from general considerations. Note that this is 



not the case in the approach of [ P4| , [351 where the results obtained with the 



Drude model do not join the results of the plasma model with decreasing 7. 
The reason is that, different definitions of the zero- frequency term of Lifshitz 
formula are used in ^ and in the present paper (see Sec. II). 



As to the papers |37|, there is the smooth transition between the 



results obtained in Drude and plasma models. However, as was indicated in 
Figs. 0,0, for a > 4/im the results in both models coincide with one another 
and for separations a > 6 fim coincide also with the case of ideal metal. 

Now let us derive some analytic results for the configuration of two plane 
parallel plates using both the plasma and Drude models. In the plasma model 
it is possible to obtain the perturbation expansion of Eq. (|231) in terms of a 
small parameter (5o/a, where 60 = c/up is the effective penetration depth of 
the electromagnetic zero-point oscillations into the metal. For this purpose 
it is useful to introduce the new variable v = ^/y instead of ^ and to rewrite 
Eq. ( |23|) in the form 



Fssia) 



he 



1^2 4^ y^dy dv COS {lvyt)fss{v,y). (52) 
iOTT a 

Expanding the quantity fss defined in ([2ID| ) up to the first order in powers of 
60/ a one obtains 



fssiv,y) = 
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ey - I 



1 + .^)^. 



(53) 
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Substituting Eq. ( |53| ) into Eq. ([5^ ) we come after some transformations to 
the Casimir force including the effect of both the nonzero temperature and 
finite conductivity 

1 TT^ cosh(7rt/) 
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1 + 



5n °° 



6 a a 



1^1 



1 



tH^ tl sinh^TTtl) 
2cosh^{7Ttl) + 1 2cosh(7rt/) 
sinh'^(7rt/) 
coth(7rt/) 



(54) 



TTtl sinh'^(7rt/) 



where Fgg(a) = F^g{a,0) = — 7r^fic/(240a'^) is the zero-temperature Casimir 
force between ideal metals. The first summation in (|5^ ) is exactly the tem- 
perature correction in the case of ideal metals (see, e.g., []T^). The second 
summation takes into account the effect of finite conductivity combined with 
nonzero temperature. 

In the limit of low temperatures T <C Tg// one has from ( |5^) neglecting 
terms exponentially small in 27TT(.ff/T 

4 



F.J a) 
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(55) 



For (5o = (perfect conductor) Eq. ( |55| ) turns into the well known result [ [TB 



demonstrating that the first nonzero temperature correction is of the fourth 
power in T/T^jf. For T = the first-order finite conductivity correction 



to the Casimir force ^n\, [TBI reproduced from (|55|). Note that the first 



correction of mixing finite conductivity and finite temperature is of order 
[T/Tgff)^. More significantly, note that there are no temperature corrections 
of order (T/Tf.ff)'^ with k < 4 in the higher order conductivity correction 
terms {So/aY from the second up to the sixth order . 



In the Drude model the analytical results can be obtained in the high- 
temperature limit T ^ Teff. It is easily seen that at high temperatures 
(large separations) only one term of Eq. (|3^) with / = contributes the 
result, the other terms with / > 1 being exponentially small in the parameter 
27TT/Teff. The situation here is exactly the same as for the ideal metal [[TBI 
As a consequence, Eq. ([351) can be rewritten in the following form 



FJa) = 
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(56) 



It is seen from Eq. ([55[) that in the high-temperature limit only the perpen- 
dicular reflection coefficient gives rise to the finite conductivity corrections 
to the Casimir force (the same is valid in plasma model also). 



After the straightforward calculations up to the first orders in small pa- 
rameters (5o/a, /ojp (see the details in Appendix II) the result 



F,,(a) = FO(a,r) 



1-3 TT^^iW 



(57) 



is obtained, where Ii is a function slowly depending on the effective relaxation 
parameter (space separation) and given by 

h{l) = Jdy T-^^- (58) 

At high temperatures, which are considered here, the Casimir force acting 
between the ideal metals is 

f.>,r) = -|^C(3), (59) 

where C(3) ~ 1.202 is the Riemann zeta function. 

With the above value of relaxation frequency for Al one has /i 1.16 at 
a = 5 fim and Ii « 0.99 at a = 10 fim. In Appendix II (Fig. g) the functional 
dependence of Ii on 7 is plotted. As is seen from (|57| ), the high-temperature 
Casimir force depends on both plasma frequency and relaxation frequency. 
For all a > 6 fim the results obtained by the asymptotic Eq. ( |57|) coincide 
with the above results of numerical computations (see Figs. and 0). The 
characteristic size of the conductivity correction at large separations can be 
estimated from the following example. At a = 10 fim the finite conductivity 
correction obtained from ([57]) is 6c{F^) « 0.89% and with 7 = it holds 
Sc{FPg) « 0.47% in perfect agreement with Figs. |2|, 0. The smooth joining 
of the results obtained using the Drude model with those using the plasma 
model when 7 ^ is evident. 

Note that in the framework of the approach used in [jS^, ^ the high tem- 
perature Casimir force between real metals is given by Fss{a) = Fgg{a,T) /2 
(see Eq. (|59| )), i.e. two times smaller than for ideal metal, regardless of 
the conductivity properties of a real metal used. As to the approach used 
in []3^, |37[], the asymptotic behavior at high temperatures coincides with 



Eq. ( |5^) , i.e. it is the same as the case for ideal metals in both plasma and 
Drude models. Once again, the real properties of a metal do not influence 
the result. 



V COMPUTATIONS OF THE CASIMIR FORCE BE- 
TWEEN A PLATE AND A SPHERICAL LENS 
MADE OF REAL METAL 

The configuration of a spherical lens (or a sphere) above a semispace (plate) 
was found to be the most suitable for the precision measurements of the 
Casimir force [j7|, g, ||, (T^, [TT|, |T2j . In these experiments the finite conductivity 
corrections have been demonstrated [g, 0, [T^, HI] the sensitivity for 
the detection of the temperature corrections is close to being achieved [j7[] . For 
this reason the combined effect of both corrections is of extreme interest. Here 
the computations of the Casimir force for the configuration of a lens above 
a plate are performed using the Drude and plasma models with parameters 
of Eq. (191). The results obtained by our Eq. (|S|) are compared with the 
computations of other authors. 

We start with the relative temperature correction 



^ Ffi{a,To)-Ffi {a,0) 



where all the notations and parameters are the same as in Eq. (^) and 
only the configuration is different. The results of numerical computations by 
Eqs. (|5|), ([5D1) in the case of Drude model are presented in Fig. |^ by the solid 



curve 1. The solid curve 2 represents the temperature correction in the case 
of a dielectric plate and a lens (see the next section). By the dashed curve the 
results obtained by the approach of |35| are shown. The curve 1 increases 
monotonically in perfect analogy with Fig. However, the dashed curve 
represents the negative temperature correction at separations a < 4.1 fim and 
changes sign for larger separations. This behavior corresponds to the case 



where large by modulus corrections linear in temperature to the ideal zero- 
temperature Casimir force are present. According to such corrections 
are in contradiction with the experimental data of [|7|. On the basis of our 
Eq. (^) such corrections do not arise. 

In Table 2 the values of the relative temperature corrections are presented 
at different distances between a lens and a plate computed using the Drude 
model (present paper, column 2) and the plasma model (present paper, col- 
umn 3), the ideal metal (column 4), computed on the basis of [j35j using the 
Drude model (column 5), and on the basis of ^ using Drude and plasma 
models (columns 6 and 7 respectively). 

From Table 2 (columns 2-4) it follows that at smallest separations the 
temperature correction computed using the Drude model is about two orders 
of magnitude larger than with the plasma model and for the ideal metal. At 
large separations the predictions of both models are very close to each other 
and to the results obtained with the ideal metal (column 4). The negative and 
large magnitude of the corrections at small separations results of column 5, 
computed on the basis of correspond to linear temperature corrections. 
The results of columns 6, 7, although positive, also correspond to the presence 
of linear temperature corrections at small separations. 

The relative finite conductivity correction at a temperature Tq can be 
described by 

^Ff,) = (61) 



In Fig. this quantity is plotted as a function of separation in the case of the 
Drude model at Tq = 300 K (solid curve is our result computed by Eq. (|48|), 
the long-dashed curve is computed on the basis of ^). The dependence 
of the conductivity correction on separation at Tq = is shown by the short- 
dashed curve. In Fig. the analogous results computed using the plasma 
model are presented. As with the case of two plane plates (Sec. IV), Drude 
model leads to larger finite conductivity corrections than the plasma model. 
At separations larger than 4 /im the conductivity correction computed using 



35|, |37] turns into zero with both Drude and plasma models (the long-dashed 



curves in Figs. 0, |7|). Once more, this property is determined by the artificial 
modification of the zeroth term of Lifshitz formula used in |37| and does 
not depend on the particular characteristics of a real metal. 

In the same way as for two plane plates our results for a lens above a plate, 
computed using the Drude model, join smoothly when 7^0 with the results 
computed with the plasma model. This is not the case in the approach of 



We come now to the perturbative analytical results which can be obtained 
for the configuration of a lens above a plate. Here, the plasma model can be 
used. Introducing the new variable v = ^/y instead of ^ one obtains 
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(62) 



The expansion of fsi up to first order in the small parameter d^/a is 
fsi{v,y) = 2\xi{l-e-y)+2- 



ey -I 

Substitution of Eq. (|63| ) into Eq. (^) leads to result 
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where ^^(a) = FS(a,0) = - tt^ hcR / {360 a^). 

For the case of low temperatures T <C Teff Eq. (^ leads to [j^^ 

3 / rr. \ 4 



F,i{a) « FO(a) 1 + 



45C(3) / T 



T 



eff, 



eff, 



60 _ 45C(3) 
'a 27r3 T, 



+ 



T 



eff, 



eff. 



(65) 



The corrections to (^5|) are exponentially small in the parameter 2TTT(,ff/T. 
For ideal metal = and Eq. ( |65|) coincides with the known result [[7j. 
It is seen that for the perfectly conducting lens and plate the first nonzero 
temperature correction is of the third order in T/Te//. For T = the first 
order finite conductivity correction to the Casimir force is reproduced. In 
analogy with two plane plates the perturbation orders {60 /aY with 2 < i < 6 
do not contain temperature corrections of orders (T/Teff)^ and [T/TpjfY or 
smaller ones [|33|] . 

Now we consider the analytical results which can be obtained with the 
Drude model at high temperature T ^ ^e//- Here only the zeroth term of 
Lifshitz formula contributes the result. Then Eq. can be represented as 



Fsi{a) = 



UbTR 
8a2 



00 



00 



ydy ln(l — e' 



ydy In [l - rl{y,y)e 



(66) 



10 



After some transformations (see Appendix II) one arrives at the result 



F4a) = F^i{a,T) 
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where I2 is defined by 



(68) 



The high-temperature Casimir force acting between the lens and the plate 
made of ideal metal is 



(69) 



For example, I2 ~ 0.519 at a = 5 iim and Ii « 0.434 at a = 10 iim (using 
the data of ([491) for ^4/). The dependence of I2 on 7 is plotted in Appendix 
II (Fig. 0). The asymptotic results of Eq. ( |67| ) coincide with the results of 
numerical computations for a > 5 fim. The choice of the model describing 
dielectric properties of a metal at large separations is rather important. At 
a = 10 /im the conductivity correction is 6c{F^) « 0.68% using the Drude 



model and Sc{Fgi) « 0.32% using plasma model, i.e. more than two times 
smaller. In analogy with two plates, the high-temperature Casimir force 
calculated on the basis of is two times smaller than in ( 1691 ) regardless 



of the conductivity properties of the metal. In |T7[] the high-temperature 



behavior for the real metals coincides with (^9]) obtained for the ideal metal. 
So, the actual properties of a particular metal are not reflected. 

VI TEMPERATURE CASIMIR FORCE BETWEEN 
THE DIELECTRIC TEST BODIES 

Here we briefly discuss the temperature Casimir force between dielectrics. 
The application of Lifshitz formulas ([I|) or (|37|) for the case of dielectric 
surfaces is direct. No additional prescriptions, such as the ones used above 
or their generalizations are needed to obtain the final result matching the 
general physical requirements. However, when the formalism developed for 
dielectrics is compared with that for metals, the origin of the above difficulties 
becomes clear. 

The dielectric permittivity of dielectrics can be modeled, e.g., by the 
Mahanty-Ninham relation [Mj 

= 1 + (70) 

-L I , ,2 

Here cjg ~ 2 x 10^^ Hz gives the main electronic absorption in the ultraviolet, 
£o is the static dielectric constant. At small ^ <C cjg one has £{i^) ~ eo- In 
fact, frequencies giving large contributions to the Lifshitz formulas or ( |37| ) 
in the micrometer separation range are much smaller than cUe- Because of this, 
£{i^) can be approximately replaced by Eq. Below we use in all computations 
£o ~ 7 which corresponds to the sheet of mica. 

In reality, the zeroth term of Lifshitz formula for dielectrics is also discon- 
tinuous as in the case of metals. To illustrate this statement we put / = 
in and obtain the following values of the reflection coefficients defined in 



Eq. (a) 

r?{0,M=(|^)', r^{0,fcx)=0. (71) 
These values do not depend on k±. Therefore they are preserved in the hmit 
k± 0. At the same time, if we put k± = from the very beginning one 
obtains 

r?(&0) = ri(&0) = rW=(^^)', (72) 

which is the case for real photons. These values do not depend on and 
are preserved in the limit 0. Eqs. ([n|), (|7^) together imply that both 
reflection coefficients ri(^, k±) and r2(^, k±) are discontinuous as the functions 
of two variables at a point (0,0). Remind that in the case of metals described 
by Drude model only the transverse reflection coefficient r2 was discontinuous 
(see Sec. II). Note that in the zeroth term of the Lifshitz formula for dielectrics 
both reflection coefficients (|7T]) correspond to nonphysical (virtual) photons 
with ri > r^^) and r2 < r^^\ At the same time, for metals the longitudinal 
reflection coefficient at zero frequency takes the physical value n = r*^^) = 1 in 
both plasma and Drude models. As to the perpendicular reflection coefficient 
at zero frequency in the case of metals, it corresponds to nonphysical photons 
and takes the values r2 = < = 1 in Drude model and r2 = g{k±) < 
j,{R) = 1 in the plasma model (see Eqs. (|T3|), ([T^ )). 

It is of paramount importance that the discontinuity of both reflection 
coefficients causes no physical problem in the cases of dielectrics. The point is 
that the dielectric permittivity ( |70| ) corresponds to a nondissipative medium 
(which is not the case for metals described by the Drude model). For this 
reason, the scattering problem, which underlies the Lifshitz theory (see Sec. II 
and also [^), is well defined at zero frequency through the unitarity of 



scattering matrix, furnishing the desired value of the scattering coefficient 
and thereby the free energy. The results obtained for dielectrics are physically 
consistent. They are immediately evident from Eqs. (|I]), ([1771) without use of 
any additional assumptions which are necessary in both cases of ideal and 
real metals where the scattering problem at zero frequency is not well defined. 



The results of numerical computations for the dielectric test bodies made 
of mica are shown by the solid curves 2 in Fig. [I| (Eq. (|I]), two plane plates) 
and in Fig. |^ (Eq. (^), a lens above a plate). The curve 2 is in direct 
analogy with the curve 1 in the same figure. At all separations the relative 
temperature correction for dielectrics is positive. At a = 0.1 fim it takes the 
value 6t{Fss) = 1-94 x 10^^. This is larger than for ideal metal and for real 
metal (^4/) considered in framework of the plasma model but smaller than for 
the same metal in framework of the Drude model (see Table I). At a = 1 fim 
and a = 10 /im there are 5t{F,,) = 1.99 X 10-2, respectively, 6t{F,s) = 2.50 
for dielectrics. At a > 1 /im the temperature correction for dielectrics is larger 
than for ideal or real metals at the same separation. For the configuration 
of a dielectric lens above a disk at separations a = 0.1, 1 and 10 /im the 
temperature correction is, respectively, 2.39 x 10^^, 6.94 x 10^^, and 4.25. 
The relationship of these values with those computed for the real and ideal 
metals (see Table II) is the same as in the case of two plane plates. 



VII CONCLUSIONS AND DISCUSSION 

As demonstrated above, the computation of the Casimir force between real 
metals at nonzero temperature is a complicated theoretical problem. The 
first contradictions between Lifshitz theory [g, [T^ applied to ideal metals 
and calculations based on quantum field theory were revealed in the 



sixties. They were resolved by Schwinger, DeRaad and Milton [JT^ by the use 
of special prescription modifying the zero-frequency term of Lifshitz formula. 
After the use of this prescription the results of [@, |T6j in apphcation to ideal 
metal became agreed with those of [j2^, ^ . 

A new stage in the solution of the problem has been started only recently. 
It was motivated by the increased accuracy of the Casimir force measurements 
and possible applications of the Casimir effect as a test of fundamental phys- 
ical theories and in nanotechnology. Different authors [j32], |37| 



applied Lifshitz theory to calculate the temperature Casimir force between 



real metals and obtained the diversed results. In ^ the Lifshitz formula 
was applied to real metals in its original form without any modification of 
the zero-frequency term. The obtained results, however, turned out to be 
in contradiction with experiment (see Sec. V) and with general theoretical 
requirements (negative temperature corrections at small distances and incor- 
rect asymptotic at high temperatures, see Sees. II, IV, V). In ^ the 
zeroth term of Lifshitz formula for real metals at nonzero temperature was 
modified according to the prescription of [|TBj formulated for ideal metals. The 
obtained results were compared with the Casimir force at zero temperature 
which was computed for real metals without use of any prescription. Such 
an approach leads to significant temperature corrections to the Casimir force 
at small separations (both in plasma and Drude models) which are linear in 
temperature and also to the absence of any finite conductivity corrections at 
the moderate separations. In ^ the computations of the temperature 
Casimir force in frames of plasma model were performed with the coinciding 
results (no linear in temperature corrections were found at small distances). 
These computations did not use any modification of Lifshitz formula. In 



the problems arising in frames of Drude model were also formulated and the 
way to their solution was directed. 

In the present paper we propose the new prescription modifying the ze- 
roth term of Lifshitz formula in the case of real metals described by Drude 
model (Sees. II, III). The necessity of this prescription has in-depth reasons 
connected with the failure of scattering formalism underlying Lifshitz for- 
mula in the case when the dielectric permittivity describes a medium with 
dissipation where the unitarity condition is absent. This prescription is the 
generalization of the Schwinger, DeRaad and Milton prescription [[TBj for the 
case of real metal. In the case of plasma model (which describes a nondissi- 
pative medium) it leads to exactly the same results as an unmodified Lifshitz 
formula. Because of this all the results obtained in [^, ^ preserve their 
validity. 



The Lifshitz formula with the modified zero-frequency term is given by 
Eq. ([351) (configuration of two plane plates) and by Eq. ( ^) (configuration 
of a lens above a plate). The detailed computations with the use of these 
equations were performed in Sees. IV, V. It was shown that the obtained tem- 
perature corrections are positive and offer the correct asymptotic behavior 
both at small and high temperatures (separations). The results obtained in 
frames of Drude model join smoothly with those obtained in frames of plasma 
model when the relaxation frequency goes to zero. The finite conductivity 
corrections to the Casimir force were computed at nonzero temperature in 
the separation range from 0.1 fim till 10 fim in frames of Drude and plasma 
models. Both the temperature and finite conductivity corrections calculated 
above possess the reasonable physical properties avoiding difficulties which 
|3^ , |35|j and in |37|. The perturbative analytical results at both 



arise m 



small and large separations are in agreement with the numerical computa- 
tions. The case of dielectric test bodies (Sec. VI) where the scattering problem 
is consistent at all frequencies illustrates the essence of difficulties arising in 
the case of real metals. 

In the nearest future one should expect the experimental registration of the 
temperature Casimir force. This will be the final answer in the discussion 
on the subject what is the temperature dependence of the Casimir force. 
Meanwhile it is necessary to proceed with the more detailed elaboration of 
the microscopic theory of dispersion forces based on quantum field theory at 
nonzero temperature in Matsubara formulation. There is a real possibility 
that the above phenomenological prescription, concerning the zero-frequency 
contribution to Lifshitz formula, will be rigorously derived, at least as a good 
approximation, in terms of the scattering theory in dissipative media. The 
final solution of this problem seems to be of large importance taking into 
account the prospective role of dispersion forces in both fundamental and 
applied science. 
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APPENDIX I 



In this Appendix Eqs. (^), ( WA ) are derived starting from Eq. (P5|). Using 
(|2I]|), Eq. (|2^) can be rewritten in the form 



Fss{a) = 



he 



y dy 



1, difUiy) 



+ 



2t 



M\y.y)-Jd^- 





:df^l\ly) 



7T 

1 



9^ 



^yfss{y,y) 



1 

+ 2 



ifssii, y) 



dUss{ly) 



:dfss{iy) , (t^ 



' 



(I.l) 



The functions and df^l^d^ are continuous at all points for both 
plasma and Drude models. Due to this it holds 



d^- 



f^l\y.y)-fil\(i,y)- 



(1.2) 



Note that the same equality is not valid for /j^^ because it is discontinuous at 
zero frequency. Substituting (|TT^) into ( pT|) and performing the cancellations 



one obtains 



Tie 



/2/'rf2/{|[/i]'(0,2/)+/i?(y,!/)] 



2t 



(1.3) 



Now let us consider two last contributions to containing the function 
A{z). With account of Eq. ([1^) r = 2TiT/Teff = 27r/t According to defini- 
tion of the integer portion function 

A =A(^=k-iikT<y<{k+ l)r (1.4) 

for A; = 0, 1, 2, ... . Using ([T^ ) the integral from (JOj) can be represented as 



OO 



y dyfss{y,y)A 

3t 



2r 
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{l + l)T 
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CXD CX3 CXD 
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OO ^ 

= E / y^dyfss{y,y). 



(1.5) 



The second integral from ( p3|) containing the integer portion function is a bit 
more complicated. Using the definition of this function, it can be represented 
in the following form 



y dy \ di 



:dUly)^(ti 



2t 




3t 



27T 



= y dy d^ 



:dfss{i, y) 



+ y dy 



2t 



2r 



+ 



{1+1)t 

+ / y'^dy 
I 9i 



9? 
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Now we calculate all integrals with respect to ^ according to 



= fss{b,y) - fss{a,y). 



(1.7) 



For all a 7^ 0, as it is in ( pB|) , Eq. (fTTI ) is valid for both polarizations (i.e. for 
both and /j^-') because the quantities under consideration are continuous. 
The result is 



oo V 
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3r 



dfssily) , (ti 



2t 



A 



27T 



y^ dy [fss{y,y) - fssir^y)] 
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2t 
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+ / y^dy [Ifssiy, y) - fss{r, y) - fssC^r, y) fss{lT, 2/)] H . 



Combining the terms with identical arguments and using (fT5|), we obtain 



jv' dy I di^Ii^A j = - t Jy' dy Mir, y) 
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+ E / y^ dy j,,{y,y). 



(1.9) 



Now we substitute ( [TBI ) and ( [T^ ) into ( P~3| ) with the result 

■fc oo 

he 



327raH 



y dy 











\ Ms\i,y) 
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fil\o,y) + fSKy.y) 

Tic °° °? 

y'dyUlr,y). (I.IO) 



With account of t = T^ff/T = nc/[2akBT) and W = ii Eq. ( pTq ) coincides 
with Eqs. ([2^), (^). 



APPENDIX II 



In this Appendix Eqs. (WA), ( W^ are obtained and the integral quantities 
Inil) {i = 1, 2) are computed. The expression under the second integral in 
the right-hand side of Eq. ( |5^ ) is approximately equal to 



y 



y 



6^ 



'^—y'^^Jy{y + i) 



fii.r 



1 ey -1 a" ' " [ey - 1)2' 

where the terms up to the first order in small parameter boja are preserved. 
Notice that the first contribution here is the same as under the first integral 
in (^). They together produce the high-temperature asymptotic for the 
case of ideal metal. The second contribution in takes the effects of 

finite conductivity into account. This second contribution can be identically 
represented as 



1)^ 



a 



ey - 1 



+ 7 



y^^/y 



\/^T7 + v^(e^-l)' 

(11.2 



Substituting with account of 26o/a = A/cUp and ^/up 

into one obtains 



F.J a) = 



167ra^ 



4C(3)-12C(3)^- 4^/1(7) 



a 



UJr 



(II.3) 



where /i(7) was defined in ([58|) . This equation coincides with Eq. ([57|) be- 
cause of (^). 

In Fig. g the dependence of I\ on 7 is plotted (curve 1). The values of 
I\ decrease from 1.3844 at 7 = 1 until 0.8782 at 7 = 10. With account of 
definition of Sec. II 7 = 2a'^ jc the dependence of 1\ on 7 may be recalculated 
into the dependence of I\ on a at fixed 7. To take an example, with the above 
value of 7 = 9.6 x lO^^rad/s 7 changes in the interval 1.92 < 7 < 6.4 when 
the separation distance belongs to the interval 3 /im < a < 10 \mi. 

We now turn to the derivation of Eq. (|67|). The expression under the 
second integral in the right-hand side of Eq. (|66|) calculated up to the first 



order in d^/a is 



.5, 



^ In 1 - r|(7/, y)e y\ ^ y ln(l - e y) + 2—^y{y + 7) — 



y 



1 



(II.4) 



Once more the first contribution here is the same as under the first integral 
in (^). Together they produce the high-temperature Casimir force between 
a lens and a plate made of ideal metal. The second contribution in ( [11 .41) is 
responsible for the finite conductivity correction. This second contribution 
can be rewritten in the form 



2— \/?/(2/ + 7) — 
a ey 



y 



6r 



y 



1 a \ey -1 

Substituting JTT^) and ( pT5| ) into Eq 
one obtains 

ksTR 
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(II.5) 



Fsi{a) 



8a2 



\/y + i + y/y - 1, 

and performing the integration 



2C(3) + 4C(3)^ + 4^/2(7) 



a 



where hil) was defined in ( |58|) . Eq. ( |[l.tj|) coincides with Eq. (f57|) if to take 
account of (|5^). 

In Fig. g the dependence of I2 on 7 is shown by the curve 2. The values 
of I2 decrease from 0.6455 at 7 = 1 until 0.3787 at 7 = 10. 
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Table 1: The relative temperature correction to the Casimir force between two plates in 
dependence of separation for different models of a metal conductivity 





Present 


paper 




Ideal metal 


Approach of 


llo rll 

i§ 


Approach of |6|, |3 


a (/im) 














X 1 I7<D\ 








U.l 


5.16 X 10" 


Q 
O 


6.60 X 10 


-6 


1.57 X 10- 


7 


—1 .1 1 X 10 


-3 


2.18 X 10"^ 


1.62 X 10 


0.3 


9.24 X 10- 


-3 


5.48 X 10- 


-5 


1.27 X 10- 


-5 


-3.49 X 10 


-2 


2.54 X 10-2 


1.53 X 10 


0.5 


1.08 X 10- 


-2 


2.11 X 10- 


-5 


9.82 X 10- 


-5 


-6.43 X 10 


-2 


2.70 X 10-2 


1.52 X 10 


0.7 


1.18 X 10- 


-2 


6.05 X 10- 


-4 


3.77 X 10- 


-4 


-9.41 X 10 


-2 


2.83 X 10-2 


1.54 X 10 


1 


1.37 X 10- 


-2 


2.06 X 10- 


-3 


1.57 X 10- 


-3 


-0.138 




3.06 X 10-2 


1.68 X 10 


3 


0.136 




0.123 




0.117 




-0.324 




0.156 


0.137 


5 


0.580 




0.563 




0. 553 




-0.185 




0.602 


0.577 


7 


1.17 




1.15 




1.14 




9.75 X 10" 


2 


1.19 


1.16 


10 


2.08 




2.06 




2.05 




0.556 




2.11 


2.07 



Table 2: The relative temperature correction to the Casimir force between a lens and a plate 
in dependence of separation for different models of a metal conductivity 
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Figure 1: Relative temperature correction to the Casimir force between two plates in 
dependence of separation. Curve 1 corresponds to Drude model (our computation), the 
dashed curve is obtained in Drude model with r2(0, k±) = 0, and curve 2 is for the dielectric 
plates. 
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Figure 2: Relative finite conductivity correction to the Casimir force between two plates 
in dependence of separation in Drude model. Solid curve represents our computations at 
T — 300 K, long-dashed curve is obtained under the supposition ri_2(0, k±) = 1 at T = 300 K, 
short-dashed curve is for T = 0. 
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Figure 3: Relative finite conductivity correction to the Casimir force between two plates 
in dependence of separation in plasma model. Solid curve represents our computations at 
T — 300 K, long-dashed curve is obtained under the supposition ri_2(0, k±) ~ 1 atT — 300 K, 
short-dashed curve is for T = 0. 




Figure 4: Relative temperature correction to the Casimir force between two plates in 
dependence of relaxation frequency at a separation a — 2 /im and T — 300 K (solid curve). 
The dashed line is for plasma model. 




Figure 5: Relative temperature correction to the Casimir force between a plate and a lens 
in dependence of separation. Curve 1 corresponds to Drude model (our computation), the 
dashed curve is obtained in Drude model with r2(0, k±) = 0, and curve 2 is for the dielectric 
test bodies. 
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Figure 6: Relative finite conductivity correction to the Casimir force between a plate and a 
lens in dependence of separation in Drude model. Solid curve represents our computations at 
T — 300 K, long-dashed curve is obtained under the supposition ri_2(0, k±) ~ 1 atT — 300 K, 
short-dashed curve is for T = 0. 
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Figure 7: Relative finite conductivity correction to the Casimir force between a plate and a 
lens in dependence of separation in plasma model. Solid curve represents our computations at 
T — 300 K, long-dashed curve is obtained under the supposition ri_2(0, k±) ~ 1 atT — 300 K, 
short-dashed curve is for T = 0. 




Figure 8: Dependence of the coefficient integrals in Eqs. (57) (curve 1) and (67) (curve 2) 
on the dimensionless relaxation frequency. 



